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Abstract. We present the first diffraction-limited images of 
the mass-loss envelope of the red supergiant star VY CMa. 
The two-dimensional optical and NIR images were recon- 
structed from 3.6 m telescope speckle data using bispectrum 
speckle interferometry. At the wavelengths ~ 0.8 /im (RG780 
filter), 1.28 /im, and 2.17 /im the diffraction-limited resolu- 
tions of 46 mas, 73 mas, and 124 mas were achieved. All im- 
ages clearly show that the circumstellar envelope of VY CMa 
is non-spherical. The RG780, 1.28 /im, and 2.17/itn FWHM 
GauB fit diameters are 67 mas x 83 mas, 80 mas x 116 mas 
and 1 38 mas x 205 mas, respectively, or 100AU x 125AU, 
120 AU x 174 AU and 207 AU x 308 AU (for a distance of 
1500pc). We discuss several interpretations for the asymmet- 
ric morphology. Combining recent results about the angular 
momentum evolution of red supergiants and their pulsational 
properties, we suggest that V Y CMa is an immediate progeni- 
tor of IRC + 10 420, a post red supergiant during its transforma- 
tion into a Wolf-Rayet star. 
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1. Introduction 

VYCMa (HD 58061, SAO 173591) is one of the most lumi- 
nous red supergiants in the Galaxy with L ~ 4 10 5 Lq (Sopka 
et al. 1995, Jura & Kleinmann 1990) and is, therefore, an ideal 
candidate for the study of the progenitor phases of a supernova. 
Its distance is about 1500 pc, its spectral type is M5 Iae, and it 
is variable with a period of ^2200 days (Jura & Kleinmann 



et al. 



199C; Malyuto et al. 11997 ; Danchi et al 



1994; Knapp & Mor- 



1985| ; Imai et al. |1997J). Herbig ( |1972| ) found that VY CMa 
is embedded in an optical nebula with a size of 8" x 12" at 
650 nm. The first high-resolution observations of the du st shel l 
of VYCMa have been reported by McCar thy & Low ( |l975| ), 
McCarthy et al. ( |l977| ) and Danchi et al. ( |l994| ). VYCMa is 
a source of H2O, OH and SiO maser emission (e.g. Masheder 
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1974, Bowers et al. 1983 



Bowers et al. 
1997|, Humphreys et al. |1997|, Richards et al. 



1993 



Imai et al. 



1998). The H 2 



masers are distributed in an east-west direction, whereas OH 
masers are distributed in a north-south direction, possibly in- 
dicating a disk and a polar outflow. HST FOC images of the 
envelope of VY CMa were obtained by Kastner and Weintraub 



(1998). They show an asymmetric flux distribution in an ap- 



proximately east-west direction with a brighter core of pure 
scattered light elongated from SE to NW. 

In this Letter we present diffraction-limited ~ 0.8 /im, 
1.28 /im, and 2.17 /im bispectrum speckle interferometry ob- 
servations of the mass-loss envelope of V Y CMa. The obser- 
vations are presented in Section 2, and the non-spherical shape 
of the envelope is discussed in Section 3. Clues for VYCMa's 
evolutionary state are derived in Section 4. 

2. Observations 

The optical and NIR speckle interferograms were obtained 
with the ESO 3.6m telescope at La Silla on February 6,7 
and 8, 1996. The optical speckle interferograms were recorded 
through the edge filter RG780 (center wavelength of the fil- 
ter/image intensifier combination: ~ 0.8 /im; effective filter 
width: ~ 0.07 /im) with o ur optical speckle camera described 
by Hofmann et al. ( 1995 ). The near infrared speckle interfero- 
grams were recorded with our NICMOS 3 camera through in- 
terference filters with center wavelength/FWHM bandwidth of 
1.28 /mi/0. 012 /im and 2.17 /irn/0. 021 /im. The observational 
parameters (number of frames, exposure time per frame, pixel 
scale and seeing) are listed in Table |} Diffraction-limited im- 
ages were reconstructed from the speckle interferograms by 
the bispectrum speckle interferometry method (Weigelt 1977; 
Lohmann et al. 1983; Weigelt 1991). The power spectrum 



of VYCMa was determined with the speckle interferome- 
try method (Labeyrie 1970). The atmospheric speckle trans- 
fer functions were derived from speckle interferograms of the 
unresolved stars H42071 (RG780, 2000 frames), IRC 30098 
(1.28 /im, 600 frames) and 1 Pup (2.17 /im, 800 frames). 

Figure |l] shows contour plots and intensity cuts of the 
reconstructed bispectrum speckle interferometry images of 
VYCMa. The resolutions of the 0.8 /mi, 1.28 /im and 2.17 /im 
images are 46 mas, 70 mas and 111 mas, respectively. The 




Fig. 1. Top: Contour plots of the RG780 (left), 1.28 /im (middle) and 2.17/im filter (right) bispectrum speckle interferometry 
reconstructions of VYCMa . The contours are plotted from 15% to 100% in 10 steps. North is at the top and east to the left. 
Bottom: Intensity cuts through the centers of the reconstructed images along the major and minor axes (solid lines). The position 
angles of these axes (RG780: 153°, 1.28 /im: 176°, 2.17/im: 160°) are taken from the two-dimensional fits to the visibility 
functions (see Tab. ||) and indicated in the contour plots. The dashed curves are cuts through the reconstructed images of the 
reference stars. 




Fig. 2. Azimuthally averaged visibility functions of VYCMa and of the corresponding model flux distribution: (left) filter 
RG 780, (middle) 1 .28 /im filter and (right) 2.17 /im filter. The diamonds indicate the observations, the solid line the azimuthally 
averaged fit curve of a two-dimensional Gaussian model flux distribution plus an unresolved object, and the dashed line the 
constant caused by just the unresolved component. The RG 780 visibility function is best described by two Gaussian flux distri- 
butions. For the fit parameters see Table ^. The visibility errors are ± 0.1 up to 80% of the telescope cut-off frequency and ± 0.2 
for larger frequencies. 



Table 1. Observational parameters. 

filter number exposure pixel scale seeing 

of frames time 

RG780 2000 6lhms 7.2 mas THi 

1.28 /im 800 150ms 23.8 mas l'.'O 

2.17 /im 800 100 ms 47.6 mas l'.'5 



envelope of VYCMa is asymmetric at each of the three 
wavelengths. The object parameters were determined by two- 



dimensional model fits to the visibility functions. The mod- 
els consist of two-dimensional elliptical Gaussian flux distribu- 
tions plus an additional unresolved component. The ~ 0.8 /im 
image is best described by two Gaussians while the 1 .28 /im 
and the 2.17/im images are best described by one Gaussian 
and an additional unresolved component. The best-fit parame- 
ters are listed in Table ||. We present the azimuthally averaged 
visibility function of V Y CMa together with the corresponding 
azimuthally averaged two-dimensional fits in Fig. ||, in order 
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Table 2. VYCMa's parameters derived from the model fits to 
the visibility functions (one two-dimensional Gaussian flux dis- 
tribution plus an unresolved object for the 1 .28 /im and 2.17 /im 
observations and two two-dimensional Gaussian flux distribu- 
tions for the RG780 data). The parameters are the position 
angle of the major axis, the axes ratio (major/minor axis), 
the FWHM of major and minor axes, the azimuthally aver- 
aged FWHM diameter and the relative flux contributions of the 
Gaussian and the unresolved object. We estimate the errors of 
the position angles to ~ ±10° and those of the FWHM sizes to 
~ ±15%. 



data set 


RG780 


1.28 ^im 


2.17 ^m 


PA (°) of the major axis 


153/120 


176 


160 


Axes ratio 


1.2/1.1 


1.5 


1.5 


Major axis (mas) 


83/360 


116 


205 


Minor axis (mas) 


67/280 


80 


138 


Average diameter (mas) 


74/320 


96 


166 


rel. flux of Gaussian 


0.75/0.25 


0.91 


0.50 


rel. flux of unres. comp. 


0.00 


0.09 


0.50 



to show the wavelength-dependent relative flux contribution of 
the unresolved component (dashed line). 



3. Interpretation 

Figures 1 and 2 and Tab. [| indicate that V Y CMa consists 
of both an unresolved component and a resolved asymmet- 
ric extended component with an increasing average diameter 
at longer wavelengths. The unresolved component is proba- 
bly the star itself or an additional compact circumstellar ob- 
ject. The relative intensity of this unresolved component de- 
creases towards shorter wavelengths. The complete obscuration 
of the unresolved component close to the maximum of the stel- 
lar spectral energy distribution indicates a high optical depth of 
the circumstellar envelope, in agreement with the value from 
Le Sidaner and Le Bertre of rio Mm = 2.4. The visibility func- 
tions, the images and the best fit parameters given in Tab. || at 
the three different wavelengths can be used as additional con- 
straints in future two-dimensional radiation transfer modeling 
of the dust distribution. The resolved structures seen in Fig. 1 
belong to the circumstellar envelope of VYCMa. In fact, the 
size of the image is of the order of the dust condensation radius 
R c which has been estimated by Le Sidaner & Le Bertre within 
their spherically symmetric model to R c ~ 12i?* with i?* ~ 
4000 i?©. They note, however, that their model fit to the spec- 
tral energy distribution of VY CMa is not entirely satisfactory 
and argue that the envelope of VY CMa may be non-spherical, 
as indicated by the high level of polarisation (Rruszewski & 
Coyne 1976). Danchi et al. (1994) found R c ~ 5i?*. 

Several interpretations for the non-sphericity seen in all 
three images (Fig. 1) are possible. The position angle of the 
major axes of the approximately elliptical shapes is similar, al- 
though not identical for all three cases (153° to 176°; see Ta- 
ble 2). 



This position angle is approximately perpendicular to the 
major axis of the H2O maser distribution (Richards et al. 1998| ) 
and similar to the distribution of the OH masers (e.g. Masheder 
et al. 1974). Accordingly, we can interprete the structure of the 
circumstellar envelope of V Y CMa as a bipolar outflow in a 
north-south direction caused by an obscuring equatorial disk 
in an east-west direction. The existence of an obscuring disk is 
supported by the obscuration of the central star at optical wave- 
lengths. Such geometry was also discussed for IRC + 10 216 by 
Weigelt et al. (1998) and Haniff & Buscher (1998) and has al- 
ready been propos ed for VYCMa due to maser observations 
by Richards et al. ( |l998 l. 

Furthermore, we can not rule out that the unresolved com- 
ponent consists of an optically thick torus hiding a close binary, 
as was proposed for the Red Rectangle (e.g. Men'shchikov et 
al. 1998). This could lead to an asymmetric outflow in a north- 
south direction. 

The mass-loss mechanism of VY CMa could also be erratic 
or stochastic, similar to the clumpy pulsation and dust-driven 
mass-loss events recently detected in the prototype carbon star 
IRC +10216 (see Weigelt et al. 1998, Haniff & Buscher 1998). 
Although the reason for this anisotropy is unknown, and the 
physics of mass-loss in oxygen-rich red supergiants may dif- 
fer from those in carbon-rich stars, the common properties of 
VYCMa and IRC +10216 — both are pulsating cool lumi- 
nous stars with extended convective envelopes — could result 
in similar mass-loss features. But individual clumps are not ob- 
servable because of the larger distance of V Y CMa, which leads 
to a regular asymmetric image elongated in a north-south direc- 
tion. 

Finally, our results are also conceivable with a more regu- 
lar geometry of VY CMa's circumstellar envelope, for exam- 
ple, a disk-like envelope which appears elongated in a north- 
south direction due to the projection angle (discussed in Sec. 4 
in more detail), which was proposed earlier for VY CMa on the 
basis of optical, infrared (Herbig 1970, 1972; McCarthy 1979; 
Efstathiou & Rowan-Robinson 1990), and maser observations 
(van Blerkom 1978, Morris & Bowers 1980, Zhen-pu & Kaifu 
1984). 



4. Evolutionary status and conclusions 

With a distance of ~ 1500pc, the luminosity of VYCMa 
amounts to ~ 410 5 L Q (see Jura & Kleinmann 1990). Fig- 
ure 3 compares this luminosity to a stellar evolutionary track 
of a 40 Mq star in the HR diagram (see Langer 1991, for de- 
tails of the computational method). It constrains the initial mass 
of VYCMa to the range 30...40M Q , in agreement with e.g. 
Meynet et al. (1994), although models with rotation — which 
are not yet available for the post main sequence evolution at 
this mass — may lead to somewhat smaller values (cf. Heger et 
al. 1997). Accordingly, it is likely that VYCMa will transform 
into a Wolf-Rayet star during its further evolution. This is sup- 
ported by the very high observed mass-loss rate of V Y CMa of 
~ 10~ 4 M & yr" 1 (Jura & Kleinmann 1990). In order to obtain 
a disk-like structure (see Sec. 3), the most likely mechanisms 
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Fig. 3. Evolutionary track of a 40 M Q star of solar metallicity, 
from the zero age main sequence to the red supergiant — > Wolf- 
Rayet star transition. The luminosity of VY CMa is indicated 
by an arrow. 

involve angular momentum. There is no indication of a binary 
companion to VY CMa, either from the spectrum of VY CMa 
(Humphreys et al. 1972) or from high resolution imaging (see 
above). While this only excludes massive companions, it ap- 
pears to be viable that the axial geometry is due to the star's ro- 
tation. Direct evidence for rotation being capable of producing 
a disk-like structure around red supergiants — possibly through 
the Bjorkman-Cassinelli-mechanism of wind compression (see 
Ignace et al. 1996) — comes from bipolar structures found in 
the AGB star V Hydrae (Stanek et al. 1995), for which rapid 
rotation (v sini ~ 13 kms -1 ) has been directly inferred from 
photospheric line broadening (Barnbaum et al. 1995). 

According to Heger & Langer (1998), red supergiants dras- 
tically increase their surface rotation rate shortly before and 
during the evolution off the Hayashi line. Therefore, a disk-like 
envelope surrounding VY CMa may indicate that such a spin- 
up is currently in progress. As strong mass-loss from a convec- 
tive star acts as a spin-down mechanism (Langer 1998, Heger 
& Langer 1998), a red supergiant must previously have lost the 
major part of its envelope for the competing spin-up process to 
dominate. This strongly supports the argument that the remain- 
ing mass of VYCMa's convective envelope is in fact small, 
and that VYCMa is just about to leave the Hayashi line. It is 
also consistent with the very long observed pulsation period of 
VY CMa of about 6 yr according to the pulsational analysis of 
Heger et al. (1997), who showed that such large periods can be 
obtained in red supergiants for small envelope masses (cf. their 
Fig. 2a). 

With this interpretation, VY CMa represents the immediate 
progenitor state of IRC +10420, currently a mid A type super- 
giant evolving bluewards on human time scales on its way from 
the red supergiant stage to become a Wolf-Rayet star (Jones 
et al. 1993, Kastner & Weintraub 1995). A comparison with 
the 40 Mq model shown in Fig. 3 predicts a current mass of 
VYCMa of 15 M Q and a surface helium mass fraction of 
Y ~ 0.40. 



It is remarkable in the present context that bipolar outflows 
are observed in IRC +10420 (Oudmaijer et al. 1994, 1996, 
Humphreys et al. 1997). A disk-like structure of VY CMa's 
envelope could be the basis for such flows, which occur when a 
fast wind originating from the star in a post-red supergiant stage 
interacts with a previously formed disk, according to hydro- 
dynamic simulations of interacting wind flows (e.g., Mellema 
1997, Garcia-Segura et al. 1998). 
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